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We present a set of equations describing the cosmological gravitational wave in a gravity theory with 
quadratic order gravitational coupling terms which naturally arise in quantum correction procedures. 
It is known that the gravitational wave equation in the gravity theories with a general f(R) term 
in the action leads to a second order differential equation with the only correction factor appearing 
in the damping term. The case for a R ab R a i, term is completely different. The gravitational wave 
is described by a fourth order differential equation both in time and space. However, curiously, we 
find that the contributions to the background evolution are qualitatively the same for both terms. 
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1. Introduction: Quantum correction procedures gen- 
erate higher order curvature coupling terms in the ac- 
tion for Einstein gravity Q . Such corrections are needed 
for more realistic understanding of the early stage of the 
universe where the curvature of the universe was large 
enough for the quantum effects, ft is known that to one- 
loop order Einstein action is modified by the additional 
quadratic order curvature coupling terms Q] . In four di- 
mensional spacetime, we have two such correction terms: 
R 2 and R ab R a b terms. The first term has attracted much 
attention in the literature in the context of cosmology 
and the cosmological structure formation theory based on 
linear gravitational stability analyses. Inflationary stage 
was first implemented in R 2 gravity before it was real- 
ized based on the scalar field Q. Now, it is well known 
that through a conformal transformation of the metric 
R 2 gravity can be transformed to a minimally coupled 
scalar field ||. Thus, despite its quadratic nature, it has 
been shown that the dynamical equations describing the 
structures in the cosmological model remain to the sec- 
ond order differential equation both for the scalar type 
perturbations and for the gravitational wave ■ Such 
a symmetry does not exist in the R ab R a b term. 

In a series of papers we will investigate the roles of the 
R ab R a b term in the context of the evolution of the cos- 
mological structures. In the present work we will derive 
the equations for the background and the gravitational 
wave. The main results are Eqs. Since the cos- 

mological stability analysis for R ab R a b term is new, we 
present the necessary details for the derivation in the Ap- 
pendices. We take MTW || convention except for setting 
c = 1 = 8ttG. 

2. Gravity theory with the quadratic order quantum 
correction terms: We considered a gravity theory repre- 
sented by the following action 



S = 



~(R+ AR 2 + BR ab R ab ) 



(1) 



where L m is the matter part Lagrangian. The gravita- 
tional field equation becomes 



Rt - ±6$R - AH {1)a b - BH {2)a b = T b a , 



H (1)a b =2R a b -25 b 1 nR- 



1 



SiR - 2RRi 



6> 



(2) 



H 



(2)a 



1 



1 



R' a , - OR" - ^SiDR + UfR^Rt + 2R cd R a 



cdb' 



In four dimensional spacetime, due to the generalized 
Gauss-Bonnet theorem, contributions from R abcd R a b c d 
term can be expressed in terms of and terms. 

These modifications of Einstein gravity were first intro- 
duced by Weyl, Pauli and Eddington Q . Not only these 
terms arise from the one-loop level quantum correction, 
these are also known to make the theory renormalizable 

B 

3. The cosmological model and the gravitational wave: 
As the model describing the background universe we con- 
sider a homogeneous, isotropic and flat model. The gen- 
eral perturbation in this background model can be de- 
composed into three types: the scalar, vector and tensor 
type perturbations. To the linear order, due to the sym- 
metry in the background, these three types of perturba- 
tions decouple from each other and evolve independently. 
In this paper we will consider only the tensor type pertur- 
bation which is transverse and tracefree; the other types 
of perturbations will be considered elsewhere. The metric 
can be written as 



ds 2 = -a 2 dri 2 + a 2 (6 af3 + 2C af} ) dx a dx 13 , 



(3) 



where a(rj) is a cosmic scale factor. C a ^(x, 77) is a trans- 
verse and tracefree perturbed order tensor variable based 
on a metric 5 a p. It represents the gravitational wave with 
— = Cp a , thus having two independent compo- 
nents which indicate the two polarization states of the 
gravitational wave. C a p is invariant under the gauge 
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transformation. The inverse metric, the connection, and 
the curvatures based on the metric in Eq. ([3]) , which are 
valid to the linear order in C a p, are summarized in the 
Appendix A. The Appendix A contains all quantities we 
need for analysing Eq. (^J) under the metric in Eq. (||) 
considering to the linear order in C a p. The Appendix 
B is another presentation using t as the time variable 
where dt = adr). rj and t are the conformal time and the 
background proper time, respectively. A prime and an 
overdot indicate the time derivatives based on rj and t, 
respectively. 

4- Equations for the background and the gravitational 
wave: Using the quantities presented in the Appendix B, 
Eq. (0) leads to: 



T ° = -3H 2 - 6 (B + 3A) (2HH - H 2 + 6H 2 H 



T" = Q 



rpOL 
1 Q 5 



(4) 
(5) 



rpOt 

1 P 



-Op 



2H + 3H 2 



+2 (B + 3A) (2H - + 12HH + 9H 2 + 18H 2 H 
+ (1 + 2AR) Dp + 2ARC% 



B{ D 



P 



3HD% - 



&{h + H 2 ^ 



A 

7/2 



A 



6(H + H 2 ) C<I-AH-C$\, 



(6) 



where R and Dp are given in Eqs. ( |25| , p2| ). The back- 
ground parts of Eqs. (|g) become (T h a = f§ + 5T§): 



-rpU 



H 2 + 2{B + 3A) [2HH - H 2 + 6H 2 H j - ■-.,„. 

H + 2 (B + 3A) (h - + 3HH + 6H< 
1 



T 



3 a 



(7) 



The second equation follows from the first one; from 
T£ a = we have T ° = -H(3T ° - T£). It is remark- 
able to see that, to the background order, contributions 
from R ah R a b term are qualitatively the same as the ones 
from R 2 term. The perturbed part of Eq. (^) becomes 

Dp* + 2A [RD'p 1 + Rdf) 



H + H 2 )D%+[H + H 2 ) C% 



= 6T$, 



D^Cp' + 3HCp!~-Cl 



(8) 



Equation (||) constitutes a fourth order differential equa- 
tion for Cp (x, t) which describes the evolution of the 
gravitational wave. 

5. Case for R 2 gravity: In the R 2 gravity we let B = 0. 
Using F=l + 2AR, Eq. (g) becomes: 



H 



F R(F-l) 
~ H F + 12F 
F-HF 1 



2F 



2F 



^ a 



J_rpO 



(9) 



The second equation follows from the first one. In fact, 
by re-interpreting F as F = 1 
becomes exactly the same as Eq. 
comes 



±B)R, Ecl © 
Equation (ft) be- 



F 



(10) 



In the large scale limit without the source term we have 
a general integral form solution 



C|(x,i) = C ^x)+^(x) 



dt 

7Jf : 



(11) 



where the integration constants c^(x) and dp(x) indi- 
cate the spatial structures in the growing mode and 
decaying mode, respectively. Apparently, in the large 
scale limit, Eq. (||) also has a general solution for Cp 
which is constant in time. Pioneering studies of the cos- 
mological gravitational wave in Eintcin gravity are pre- 
sented in |llf| . Einstein gravity corresponds to a limit of 
A = = B, thus F = 1. The equations for the back- 
ground and the gravitational wave which are valid in a 
more general context of /(</>, R) (cj> is a scalar or dilaton 
field) gravity were derived in Bh see Eqs. (26,39) in Q| 
and Eq. (102) in H. Equation (10) remains valid in these 



types of more generalized gravity theories; compare with 
Eq. (26) in [Q. Cosmological studies in these types of 
generalized gravity theories can be found in [|| . 

6. Discussions: We have derived the equations de- 
scribing the evolution of the background universe and 
the gravitational wave in the gravity modified by consid- 
ering the quadratic order correction terms in curvatures. 
Such correction terms generically appear as the one-loop 
quantum correction in the investigation of the quantum 
aspects of the gravity It is known that general f(R) 
term in the action can be transformed to the Einstein 
action by a conformal transformation of the metric. The 
basic reason why we obtain a second order differential 
equation for the gravitational wave in the R 2 gravity may 
be traced to the underlying symmetry under the confor- 
mal transformation, 1^,0]. However, R ab R a i, term cannot 
be transformed to the Einstein one. 

Equation (Q) describes the evolution of the background 
universe, and Eq. (||) describes the evolution of the gravi- 
tational wave under such a background model. Although 
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the contributions to the background model are qualita- 
tively the same for R 2 and R ab R a b terms, contributions to 
the evolution of the gravitational wave are different. The 
similar contributions to the background equation are, in 
fact, an unexpected and nontrivial result. Wc hope to 
pay more attention to consequences of Eqs. in fu- 

ture investigations. 

Recently the string theories attracted much attention 
as the potential candidates for a successful quantum grav- 
ity. The generic low energy effective action of the string 
theories is represented by an action modified by a dilaton 
field. The equation for the gravitational wave in such a 
case is described by an equation which is effectively the 
same as Eq. (|To|), see j|. However, the stringy quan- 
tum correction terms are more general than our action 
in Eq. (^) ; it yields coupling terms involving various com- 
binations of the curvatures, the dilaton field, and their 
derivatives Jl2| . Similar investigations in such a more 
general gravity theories will be considered in future. 

Results for the scalar and vector type perturbations 
and applications of our results to the relevant cosmolog- 
ical situations will be presented in future publications. 

We thank Prof. S.-J. Rey for his thoughtful comments. 
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APPENDIX A: LINEAR ORDER QUANTITIES 



2—Cf - AC| 



The inverse metric (0 = rj): 

g m = -^, 5 0q = 0, g^ = \ (S a P-2C a P). (12) 



The connections: 



r[!o — —j r[5 Q — o, — + + 2— -C Q; 3, 



v a — n r a - — a q -i- n a ' 



1 fry - U /3,-y ■+■ °7,/3 °/3 7 • 

Riemann curvature: 

R bOO = = R°0ab' 

R °aO0 = ( — ) <W + Cq/3 + ~^a0 + 2 ( T" ) 

p0 /"f/ _ /T/ 

^ q/3 7 — °q 7 ,,9 u o(?,7i 

00,3 ~ 1 ~a ) 13 13 ~a 13 ' 

pa s~ia/ r~ia/ pa /la/ /nr' | Q 

21 0/37 — ~~ °/3»7> 11 /307 — °7,/3 — °/37 ' 



(13) 



R a 



0-yS 



Curvature combinations 

,\ 12 



no nfc 12 

a 4 



a \ /a 



cdO 



R cd R° 



R R 



— R cd R a cd0 , 




Covariant derivatives: 



R ° n — o [ R" R' ) i R'°m — — R' 

a A \ a 

R ' a f) = -tf (^ S P + C f ) l <"- 



OR 



~[R" + 2-R' 



(21) 



□i?° = -if<' + 2^ 



-i 2 i ^ 



OR° = = OR£, 



2— L> 



2 ( - ] D| - AD£ 



1 a' 



1 



where 



,2 l°/3 



2-C7|' - AC| 



(22) 



(23) 



and an overbar indicates a background order of the vari- 
able. 

APPENDIX B: IN TERMS OF t 

Using t as the time variable we have (H = a/a): 
i?g = 3 (i? + H 2 ) , R° a = = R%, 



Ha 



H + 3H 2 ) 8% + C% + mC% - -jC|, 



R = Q[H + 2W 



u 



-3H 1 



= G£, 



RlR b a = 12 (h 2 + 3HH 2 + 3iJ 4 



(24) 
(25) 

(26) 
(27) 



R cd R° 



cdO 



-3[H + H 2 ) (H + 3H 2 



pea dU n T?cd oa 

H H cda ~ U - H H cdOi 

R cd R a cdi3 = - (ZH 2 + 8HH 2 + 9H 4 ) 8$ 

m + 2H 2 ^ - (jH + 6i? 2 ) HCf 3 



H + 2H 2 ) ^C%, 



A 

a 5 



R' J 



-R, R' a — — R^oi 



R\ = - (H5 a p + Cp*)R, 



UR = -[R+ 3HRJ , 



(28) 

(29) 
(30) 



OR° = -3 (H- + 5HH + 2H 2 + 2H 2 H^j , 
UR° a = = aR%, 

□ = -8% (h- + 9HH + 6H 2 + 22H 2 H 
D% + 3HD" - 2H 2 D% - 



-8HHC%, 



A 



D% = C$ + 3HC%-- s qt. 
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